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Abstract. Nulling interferometry, a powerful technique for high-resolution imaging of the close neighbourhood of bright astro- 
physical objets, is currently considered for future space missions such as Darwin or the Terrestrial Planet Finder Interferometer 
(TPF-I), both aiming at Earth-like planet detection and characterization. Ground-based nulling interferometers are being stud- 
ied for both technology demonstration and scientific preparation of the Darwin/TPF-I missions through a systematic survey 
of circumstellar dust disks around nearby stars. In this paper, we investigate the influence of atmospheric turbulence on the 
performance of ground-based nulling instruments, and deduce the major design guidelines for such instruments. End-to-end 
numerical simulations allow us to estimate the performance of the main subsystems and thereby the actual sensitivity of the 
nuUer to faint exozodiacal disks. Particular attention is also given to the important question of stellar leakage calibration. This 
study is illustrated in the context of GENIE, the Ground-based European Nulling Interferometer Experiment, to be installed at 
the VLTI and working in the L' band. We estimate that this instrument will detect exozodiacal clouds as faint as about 50 times 
the Solar zodiacal cloud, thereby placing strong constraints on the acceptable targets for Darwin/TPF-I. 

Key words. Instrumentation: high angular resolution - Instrumentation: interferometers - Techniques: interferometric - 
Circumstellar matter - Planetary systems 



• '— j 1. Introduction The princip l e of nulling interferometry, first proposed 

K*N , . ... . . by iBracewelll (Il978h and general i zed by several au- 

H : The existence of extraterrestrial life in the Universe is a thors (lAnge]lll989t iLeger et alJll996t lAngel & WooHl[l997t 
. - . long-standing question of humankind. Fumre space missions fearlsson & Mennessonll2000t lAbsil et alJl2003cl) . iTto com- 
such as Darwin ( Fridlund 2000) or theTerrestrial Planet Finder ^j^^^ ^-^^^ collected by two or more telescopes in a co- 
Interferometer Peichmanetalj M | Lawson & Dooleyl axial mode, adjusting their respective phases in order to pro- 
|2005|) are currently being studied respectively by ESA and ^^^^ ^ ^^^^jjy destructive interference on the optical axis. 
NASA to search for evidence of life on Earth-hke exoEto- ^^e interferometer is characterized by its transmission map 
ets by means of nulling interferometry ^ Bracewell 1197^). r,(0, 0), displayed in Fig. |I]for a two-telescope interferometer, 
These ambitious missions ai-e particularly innovative from ^ ^^-^^^ ^ pj^^^^ ^^^^^ ^ ^^^-^^^ j^^^ ^^^^ 
the technology point of view and demand several tech- ^^^^^ (Bracewell interferometer). Its expression results from 
niques to be first tested on ground. Both ESA and NASA ^^e combination of a single telescope diffraction pattern with 
have therefore initiated studies for ground-based infrared ^j^^ ^^-^^^ ^^^^^^ produced by the interference between the 
nulhng interferometers (or "nuUers"). to be_mstalled re- beams, i.e., for a two-telescope interferometer: 
spectively at t he VLTI (iG lindema nn et alJ |2004) or Keck 

Interferometer dColavita et aL .20041) . A crvogenic nulling l2Jx{n6DI A)^ ■ il BO \ 

beam-combiner is also under develop ment for the LBT in th e T,i{0,4>)^\ ^gj^^^ j sin |^7ry cos 0j , (1) 

context of the NASA Origins Program jHerbst & Hinzl2004 . 

where 6 and (p are respectively the radial and polar angular co- 

Send offprint requests to: O. Absil ordinates with respect to the optical axis, B the interferometer 

* O. A. acknowledges the financial support of the Belgian National baseline, D the telescope diameter and A the wavelength. In the 

Fund for Scientific Research (FNRS). following study, we assume that recombination and detection 
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are both done in the pupil plane, so that no image is formed: 
the flux is integrated on the field-of-view and detected on a sin- 
gle pixel. Performing the detection in an image plane would not 
significantly change the conclusions of this study. 
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Fig. 1. Monochromatic transmission map for a 47 m Bracewell 
interferometer formed of two 8-m telescopes. The interfero- 
metric field-of-view, limited by the use of single-mode waveg- 
uides (see Sect. 13. 3> . is given by Q = jS with S the telescope 
surface. Its diameter is respectively of 137 mas and 378 mas at 
the centers of the L' and N bands. The transmission map acts as 
a "photon sieve": it shows the parts of the field that are trans- 
mitted (white stripes) and those that are blocked (dark stripes, 
including the central dark fringe) by the interference process. 
The whole transmitted flux is integrated on a single pixel. 

Besides pure technology demonstration, ground-based 
nuUers have an important scientific goal in the context of the 
Darwin/TPF-I missions, through the detection of faint circum- 
stellar dust disks ("exozodiacal dust disks") around nearby 
main sequence stars. Such dust clouds may present a severe 
limitation to Earth-like planet detection with Darwin/TPF-I 
if they are more than 20 times as dense as the Solar zodia- 
cal cloud (.Absil ^OOl) . A systematic survey of nearby main- 
sequence stars is thus mandatory to select convenient targets 
for these future space missions. Photometric surveys with in- 
frared space telescopes such as IRAS, ISO and Spitzer have 
detected cold dust around nearby stars at a level of ~ 100 times 
our Kuiper Belt, but did not allow the characterization of 
warm dust emission at a better level than ab out 500 times the 
Solar zodiacal cloud dBeichman et aklEoOsI) . Attempts to spa- 
tially resolve faint exozodiacal clouds with single dish tele- 
scopes inJhejmd4nfrM;ed_^^^ and near in- 
frared dKuchner & Brownl200ol) have not yielded better detec- 
tion limits. Even the MIRI instrument onboard the future James 
Webb Space Telescope, equipped with coronagraphic devices 
at 10.6 ;um and other longer wavelengths, does not have a good 
enough angular resolution (~ 330 mas) to detect warm dust 
within the habitable zone of Darwin/TPF candidate targets. 

Nulling interferometry, combining high dynamic range and 
high angular resolution, is particularly appropriate to caiTy out 
this systematic survey of the Darwin/TPF-I candidate targets. 



The Darwin target catalogue inclu des F, G. K and M-type main 
sequence stars up to 25 pc ( Stank ov et al. Eo05). with median 
distances to these four stellar types of respectively 20, 20, 19 
and 15 pc. Ideally, the survey should detect exozodiacal clouds 
about 20 times as dense as the Solar zodi acal cloud ("20-zodi " 
clouds hereafter). Based on the model of iKelsall et alJ (Il998l) . 
such clouds are typically 10^ times fainter than the star in the 
L' band (centered around 3.8 fim) or 10^ times fainter in the 
N band (centered around 10.5 jjm). The M band will not be 
considered in this study, because is too much aff ected by wa- 
ter va por to be useful for nulling interferometry (lYoung et alJ 
11994 . 

Besides this survey aspect, ground-based nullers have var- 
ious interesting scientific applications. Any target requiring 
both high dynamic range and high angular resolution will ben- 
efit from their capabilities, including the detection and char- 
acterization of debris disks around Vega-type stars, proto- 
planetary disks around Young Stellar Objects, high-contrast bi- 
naries, bounded brown dwarfs and hot extrasolar giant plan- 
ets. Extragalactic studies may also benefit from ground-based 
nulling through the characterization of dust tori around nearby 
AGNs. Ground-based nullers are extremely powerful at angular 
diameter measurements as well: thanks to their very high sta- 
bility, they could, for instance, detect pulsations in unresolved 
stars that still produce a non negligible stellar leakage in the 
nulled data. 

In this paper, we estimate the performance of ground-based 
infrared nulling interferometers for the detection of exozodia- 
cal dust around typical Darwin/TPF-I targets. The structure of 
the paper follows the path of the light, beginning with the as- 
trophysical sources and ending up after detection with the data 
analysis procedures. In section |3 we present the different con- 
tributors to the detected signals in ground-based nulling inter- 
ferometry. Section |3l then focuses on a particularly important 
contributor (instrumental leakage), which strongly depends on 
the turbulent processes taking place in the Earth's atmosphere. 
We will show that this stochastic contributor needs to be com- 
pensated in real time at the hardware level in order to reach the 
required sensitivity. This has important consequences for the 
design of a ground-based nulling interferometer, as discussed 
in section|3 where real-time control loops are addressed. After 
passing though the instrument and being recorded by the de- 
tector, the signal still needs to be carefully analyzed in order to 
extract the useful exozodiacal contribution from the raw data. 
Section|5]describes the associated post-processing techniques. 
In all these sections, we will take a Sun-like G2V star at 20 pc 
as a typical target for illustrative purposes. Finally, these con- 
cepts are applied in section |6l where we derive the expected 
sensitivity of a specific instrument to the dust emission for dif- 
ferent types of targets in the Darwin star catalogue. 

2. Signals and noises in ground-based nulling 
interferometry 

2.1. Exozodiacal dust disk 

The Solar zodiacal cloud, a sparse disk of 10-100 fim diame- 
ter silicate grains, is the most luminous component of the solar 
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system after the Sun. Its optical depth is only ~10"^, but its in- 
tegrated emission at 10 fj.m is about 300 times larger than the 
flux of an Earth-sized planet. Very little is known about dust 
disks around main sequence stars. We will thus assume that 
exozodiacal clouds have the same d ensity distribution a s the 
Solar zodiacal cloud, as described by Kel sall et alJ ( 1 1998 ). but 
with a different scaling factor In the following study, we will 
also assume exozodiacal disks to be smooth (no clump or wake 
in the disk) and to be seen face-on. With such assumptions, the 
azimuth of the baseline has no influence on the observations: 
only the baseline length matters, which is convenient for illus- 
trative purposes. Nevertheless, we should note the influence of 
the disk shape on the transmitted signal. 

- The disk inclination has a strong influence on the transmit- 
ted amount of exozodiacal light. If the disk is seen edge- 
on and oriented perpendicularly to the fringe pattern (i.e., 
major axis of the disk parallel to the baseline), it will pro- 
duce almost the same transmitted signal as a face-on disk. 
However, if it is oriented in the same direction as the fringe 
pattern, a large part of the disk emission will be cancelled 
out by the central dark fringe. The azimuth of the base- 
line thus has a large influence on the detection of inclined 
disks, which is an advantage: if one lets the projected base- 
line evolve as an effect of the Earth diurnal rotation, the 
exozodiacal disk signal will be modulated as the baseline 
orientation changes, making the identification of the disk 
signal easier. The effect will be all the la rger when the disk 
inclination is close to edge-on JAbsil et alJ2003bh . 

- Inhomogeneities in the disk could lead to unexpected mod- 
ulation as the baseline changes. In fact, the exozodiacal dust 
emission is expected to be smooth (less than 1% random 
variations) except for rings and wakes due to gravitational 
trapping b y planets or bands du e to recent asteroid or comet 
collisions jPermott et al.ll998l) . The detection of such fea- 
tures would be a direct signature of the presence of planets 
or planetesimals in the disk. This would however require 
a good signal-to-noise ratio since these structures will be 
much fainter than the integrated disk emission. 

2.2. Geometric stellar leakage 

Even when the star is perfectly centered on the optical axis, a 
part of the stellar light still leaks through the transmission map 
due to the finite extent of the stellar disk, an effect known as ge- 
ometric stellar leakage. We define the nulling ratio N as the ra- 
tio between the transmitted stellar flux and the total stellar flux 
collected by the two telescopes. Assuming that the stellar an- 
gular radius 0» is small as compared to the fringe spacing A/B, 
and computing the transmitted flux as a two-dimension inte- 
gral of the transmission map on the stellar disk, one gets the 
following expression for the nulling ratio: 




assuming a uniform brightness across the stellar disk. The in- 
verse of this quantity is called the rejection rate p, which is 
the fundamental figure of merit for a nulling interferometer. 



Typical values of the rejection rate for a Bracewell interferom- 
eter observing a Sun-like star at 20 pc are given in Tabled 



Baseline 


14.4m 


47m 


85m 


L' band (3.8 pm) 


22197 


2084 


637 


N band (10.5 /im) 


169573 


15897 


4863 



Table 1. Rejection rate for a Bracewell interferometer observ- 
ing a Sun-like star at a distance of 20 pc in the L' and N bands, 
for typical interferometric baselines (14.4 m - LBT, 47 m = 
VLTI, 85 m = Keck). The star is assumed to be at zenith. 



The rejection rate is quite sensitive to limb darkening, be- 
cause it is precisely the outer parts of the stellar disk that mostly 
contribute to stellar leakage. Assuming a simple linear limb- 
darkening law B^ijj) = B^i(l)[l - M^(l - yu)] for the stellar sur- 
face brightness B^, with fi the cosine of the angle between the 
normal to the surface and the line of sight, the nulling ratio 
would be the following: 




showing a linear dependance with respect to the limb- 
darkening parameter u^. For a typical K-band limb darkening 
mk = 0.26 applicable to a Sun-like star (Claret et al. 1995), the 
rejection rate would increase by a substantial 14% assuming 
the same physical diameter in the uniform and limb-darkened 
cases. Other potential contributors to the actual amount of stel- 
lar leakage are the asymmetry of the stellar disk, through stel- 
lar oblateness and/or spots. However, these are only second or- 
der effects for late-type main-sequence stars, which are gener- 
ally close to spherical symmetry and are expected to have only 
small and scarce spots. These effects will be neglected in the 
discussion. 

Due to the limited rejection rate, geometric leakage will of- 
ten exceed the flux from the exozodiacal disk at the destructive 
output of the nuller, especially in the L' band. Fortunately, geo- 
metric leakage is mostly deterministic, allowing the prediction 
of its contribution. This calibration procedure is discussed in 
Sect.ES 

2.3. Instrumental stellar leakage 

The expression (|3j of the rejection rate is valid only for a per- 
fect Bracewell interferometer. In practice, the rejection rate is 
degraded by atmospheric turbulence and various instrumental 
effects causing imperfect co-phasing of t he light beams, inten- 
sity mismatches and polarization errors ("oil ivier '199^. This 
contribution, called instrumental leakage, adds to the geomet- 
ric leakage at the destructive output of the interferometer It is 
dominated by non-linear, second order error terms and does not 
depend on the stellar diameter to the first order (Lav 2004). 

Instrumental leakage has two effects on the performance 
of the nuller: first, it introduces a bias, the mean contribution 
of instrumental leakage, and second, it introduces an additional 
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stochastic noise through its fluctuations. While the second con- 
tribution can be reduced by increasing the observation time (see 
Sect. 14. 4L the first one does not improve with time and limits 
the global performance of the nuller Unlike geometric leak- 
age, an analytical expression of instrumental leakage cannot be 
obtained: it depends on the particular shape and amplitude of 
the power spectral densities of various atmospheric effects and 
instrumental errors such as piston, dispersion, wavefront er- 
rors, polarization errors, etc ( Lav 2004). Its reduction can thus 
only be done through real-time control of these errors. In prac- 
tice, instrumental leakage can be at least partially calibrated 
by observing "calibrator stars" with well-known diameters, a 
common method in stellar interferometry. This calibration, dis- 
cussed in Sect. 15.31 will help remove the bias associated with 
the mean instrumental leakage and thereby relax the require- 
ments on the performance of the real-time control loops. 

2.4. The sky background 

The thermal emission of the sky, which peaks in the mid- 
infrared, is obtained in good approximation by multiplying the 
Planck blackbody function by the wavelength-dependent sky 
emissivity. The infrared sky brightness is plotted in Fig.|2for 
an atmosphere at a mean temperature of 284 K, which is typi- 
cal for Cerro Paranal. The mean sky brightness in the relevant 
atmospheric wavebands is given in Table0] 




2 4 6 8 10 12 14 

Wavelength [^m] 



Fig. 2. Infrared sky brightness for typical Paranal atmospheric 
conditions. 



The fluctuation of the infrared background radiation pro- 
duced by the atmosphere (or sky noise) has not been wefl mod- 
elled to date. Sky noise measurements have been carried out in 
the N band at Paranal w ith the VLTI mid-infrared instrument 
MIDI jAbsilet allbool . while L'-band measurements have 
been carried out at Siding Springs, Australia (I Allen & BartonI 
Il98ll) . showing that, at low frequencies, sky noise dominates 
over the shot noise associated with the mean background emis- 
sion. Typical logarithmic slopes ranging from -1 to -2 have 
been recorded for the sky noise power spectral density (PSD). 



2.5. The instrumental background 

In order to conduct the light beams from the telescopes to the 
instrument, a large number of relay optics are used. All of them 
emit at infrared wavelengths. The wavelength-dependent emis- 
sivity of the optical train is theoretically equal to the comple- 
mentary of its transmission. In order to evaluate this contribu- 
tion, we have to make some assumptions about the interfero- 
metric facility and nulling instrument. We will take the case 
of the VLTI, which has a transmission ranging between 30% 
and 40% in the near- and mid-infrared', and we will assume 
that the nulling instrument is at ambient temperature, except 
for the detection unit (spectrograph and detector). The instru- 
mental background is given in Tabled] separated into "VLTI 
brightness" for the interferometric infrastructure and "GENIE 
brightness" for the nulling instrument, assuming an emissivity 
of 30% for the latter. A full cryogenic instrument would reduce 
the background emission, and thus improve the global sensitiv- 
ity. However, even at ambient temperature, the full potential of 
the instrument can still be achieved by increasing the integra- 
tion time, thereby reducing the contribution of the shot noise 
associated with the background emission. Fluctuations of the 
instrumental background are also expected to happen (e.g. due 
to beam wandering, mirror vibrations, gain fluctuations, etc.), 
but at a much lower level than sky noise. They will be neglected 
in the following discussion. 

2.6. Detection criterion 

In addition to the classical sources of noise such as shot noise, 
read-out noise and background noise, whose effect can be re- 
duced by increasing the integration time, nulling interferometry 
is faced with two types of biases that must be removed to ensure 
a secure detection of the exozodiacal cloud. These biases are 
associated with the mean contribution of geometric and instru- 
mental stellar leakage. In order to detect circumstellar features 
as faint as 10"^ (N band) to 10""* (L' band) of the steflar flux 
with a signal-to-noise ratio of 5, one should reduce the con- 
tributions of both geometric and instrumental stellar leakage 
down to 10"^ (A^ band) or 10"^ (L' band) of the initial stellar 
flux, so that the total stellar leakage does not exceed about one 
fifth of the exozodiacal disk signal. Real-time control systems 
at hardware level and calibration techniques at software level 
will be used to reach this level of performance. The require- 
ments on phase and intensity control will be particularly tight 
for short operating wavelengths. In the next section, the discus- 
sion of atmospheric turbulence is thus mainly illustrated in the 
L' band where its effects are the largest. 



3. Influence of atmospheric turbulence 

In order to estimate the performance of a nulling interferometer 
in the presence of atmospheric turbulence, we have developed 



' Due to the additional emission of dust particles located on the 
optics, the measured emissivity of the VLTI optical train is in fact 
closer to 100% than the expected 60-70%. 
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a software simulator called GENIEsim JAbsil et alJ 12003 j) . 
This semi-analytical simulator computes in real-time the trans- 
mission map of the interferometer, taking into account phase 
and intensity fluctuations associated with atmospheric turbu- 
lence. It inhe rently takes into account the bilinear error terms 
discussed by La\ that are supposed to be the largest 

contributors to the instrumental nulling. In the following para- 
graphs, we describe the various sources of phase, intensity and 
polarization errors that have been implemented in the simulator 
and discuss their effect on instrumental nulling. 

3.1. The piston effect 

The Oth-order term of atmospheric turbulence, the spatially 
averaged phase perturbation on the pupil, is known as the 
piston effect. Optical path diff'erence (OPD) fluctuations are 
mainly due to differential piston between the two apertures 
of the interferometer. Under the Taylor hypothesis of frozen 
turbulence, where a frozen phase scr een is blown at the wind 
speed above the telescopes (Roddief Il98ll) . the classical de- 
scripti on of atmospheric turbulence proposed by iKolmogorovl 
( Il94ll) states that the temporal PSD of phase fluctuations fol- 
lows a -8/3 power-law. In practice, three physical phenom- 
ena in duce a departure f rom the theoretical -8/3 logarithmic 
slope dConan et al ]ll995h : correlations between the pupils and 
the finite outer scale of turbulence Xo both reduce the fluctua- 
tions at low frequencies, while pupil averaging reduces them at 
high frequencies (Fig.|3}. The latter effect has in fact never been 
observed, so that we keep the -8/3 slope at high frequencies as 
a conservative scenario. This scenario also features additional 
instrumental contributions: 



E 



--2/3 


1 ' ' 


1 1 1 




"""-■•,-3/3 






1 . < 


'■\_-17/3 \. 



0.01 /■ 0.1 /2 1 10 100 '000 



Frequency [Hz] 

Fig. 3. Differential piston power spectrum, under the following 
conditions: Z? = 47 m, D = 8 m, v = 11 m/s and = 10 cm at 
500 nm (equivalent to 1" seeing). The dotted line represents 
the theoretical Kolmogorov piston PSD, showing the effects 
of pupil correlations (below /i ^ 0.2v/B) and pupil averag- 
ing (above /2 - 0.3v/D), but without the effect of the outer 
scale nor the instrumental contributions. The solid line is the 
representation of the piston PSD that we actually use in our 
simulations as a conservative scenario. It includes the effect 
of pupil correlations and of the outer scale of turbulence (flat- 
tening of the PSD for frequencies below v/Xo), as well as the 
contributions of bimorph and coupled piston (high-frequency 
fluctuations), but excludes pupil averaging. A best-case sce- 
nario would coiTespond to the combination of the dotted line 
for / > /2 and the solid line for f < fi- 



- Bimorph piston effect: this is the OPD induced by 
resor iances of the adap tive optics deformable mir- 
ror dVerinaud & Cassaing I [2001). This contribution has 
been evaluated for the MACAO adaptive optics at VLTI, 
and turns out to be dominant above 100 Hz despite the 
use of a "piston-free" algorithm for the control of the de- 
formable mirror. 

- Coupled piston: this is the additional OPD created by 
high-order wavefront err ors when injecting the bea ms into 
single-mode waveguides dRuilier & Cassainj200lh . 

The standard deviation of the OPD fluctuations, neglecting 
the effect of the outer scale, is given bv (.Roddier.l98ll) : 

A I B\^^^ 

(roPD=2.62— - . (4) 
27r \ro/ 

Since the Fried parameter ro is proportional to A^^'^, ctopd does 
not depend on wavelength. For a typical baseline of 47 m and 
1" seeing, the theoretical variance of piston is about 35 jum, i.e., 
about 9 fringes in the L' band. This variance, which is also the 
area under the dotted curve in Fig.|3l sets the level of the PSD. 

The effect of phase eiTors on a Bracewell nulling interfer- 
ometer is to shift the position of the fringes in the transmission 

^ GENIEsim, written in IDL, is open source and available on de- 
mand. Send request to R. den Hartog (rdhartog@rssd.esa.int). 



map: the dark fringe is not centered any more on the optical 
axis. For a small phase error e^{A) between the two beams, the 
shifting of the dark fringe produces a non-null on-axis trans- 
mission r,)(0,0) = e^iAflA- (Oflivier 1999). In order to keep 
the instrumental rejection rate above the required 10^ in the L' 
band, the phase error should be smaller than 0.006 radian (i.e., 
4 nm in the L' band), while atmospheric piston induces OPD 
eiTors of about 35 jum. A deep and stable null can thus only be 
achieved by stabilizing the dark fringe on the optical axis, by 
means of a fringe tracker (Sect. RTTl . 

3.2. Longitudinal dispersion 

In addition to the achromatic piston effect produced by the fluc- 
tuations of the air refractive index, another source of phase er- 
rors comes from the fluctuation of the water vapor column den- 
sities above the two telescopes. A useful concept to describe 
the column density fluctuations of water vapor is "Water- vapor 
Displacing Air" (WDA), introduced by i Meisner & Le Pool3 
(2003). The fluctuation of the WDA column density induces 
a chromatic OPD difference between the beams, because its 
refraction index hwda strongly depends on wavelength in the 
infrared as illustrated in Fig.0] This effect, refeiTed to as longi- 
tudinal dispersion, has the same kind of influence on the trans- 
mission map as the piston effect, except that the offset of the 
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8 10 
Wavelength [fim] 



Fig. 4. Reduced refraction index of WDA(nwDA = ^^^^^aJ., with 
c the speed of light and p the molar density) expressed in fem- 
tosecond per mole/m^, plotted at infrared wavelengths. 



dark fringe is now wavelength-dependent. This has two delete- 
rious consequences on the performance of a nulling instrument: 

- Inter-band dispersion is the phase error measured at the 
center of the science waveband assuming that piston is per- 
fectly corrected at another wavelength by means of a fringe 
tracker (this correction is generally done in the H and/or K 
band, e.g. at the Keck and VLT interferometers), 

- Intra-band dispersion refers to the differential phase error 
within the science waveband: if one manages to co-phase 
the beams at the central wavelength of the science wave- 
band, their phases will not be perfectly matched at the edges 
of the band. 

Since longitudinal dispersion is produced by the same 
turbulent atmosphere as piston, it is expected to follow the 
same theoretical PSD. This has been confirmed experimen- 
tally bv Lay, ( 199Z) . An RMS value of 1.5 mole/m^ for the 
col umn density fluctuation of w ater vapor has been measured 
by 'Meisner & Le Poolel ll2003l) at Paranal on an integration 
time of 100 sec (this value corresponds to 27 yum RMS of dif- 
ferential precipitable water vapor). This translates into typical 
inter-band and intra-band OPD errors of 550 nm and 150 nm 
RMS respectively in the L' band assuming that fringe tracking 
is performed in the H band. A control system must thus be de- 
vised to stabilize the dark fringe on the optical axis in order to 
meet the requirements on the instrumental rejection rate. 

3.3. Wavefront errors 

Propagation through the turbulent Earth atmosphere degrades 
the shape of the wavefronts, producing imperfect interference 
between the beams at recombination, and thus creates a halo of 
incoherent stellar light at the detector. Modal filtering by^meari s 
of single-mode fibers, as proposed by'Mennesson et al." (2002), 
eliminates wavefront corrugations by projection on the funda- 
mental mode of the fiber, thereby ensuring perfect matching 
of the wavefronts over a broad bandpass. After injection into a 
single-mode fiber, the shape of the initial wavefront only aff'ects 



the amount of energy coupled into the guide. The modal filter 
thus converts phase defects into inte nsity errors, which are less 
severe for a nulling interferometer dOlliviei^ 1999). It also in- 
duces a small coupled piston due to the non-null mean of z > 1 
Zernike polynom ials projected onto the Ga ussian fundamental 
mode of the fiber jRuilier & C assaing 2001). The use of single- 
mode fibers to filter the beams at the focus of the telescopes 
strongly relaxes the constraints on pointing accuracy and low 
order aberrations, and thus improve the capability to achieve 
deep rejection ratios. We will assume that modal filtering, or at 
least spatial filtering using pinholes for the N band, is used in 
our nulling interferometer 

In the following study, we will also assume that the tele- 
scopes are equipped with adaptive optics, with performance 
similar to t hat of the MACAO systems on ESO's Very Large 
Telescopes jArsenault et al...2003.) . The power sp ectral densi- 
ties of individual Zernike modes are taken from IConan et alJ 
(Il995t) . with typical high-frequency slopes of -17/3 and var- 
ious slopes at low frequencies. These PSDs are corrected by 
the adaptive optics servo system for frequencies below about 
10 Hz, thereby making their low-frequency content negligible. 
The case of tip-tilt is illustrated in Fig. |5l including the con - 
tribution of high order Zernike modes qten Brummela ail 1 995t) . 
The values for the RMS tip-tilt and Strehl are given in Table|2] 

TIp-+Ilt power spectrum 




0.1 1.0 10.0 100.0 1000.0 

Frequency [Hz] 

Fig. 5. PSD for tip-tilt before (dotted line) and after (solid line) 
close-loop control with MACAO. The peaks at high frequen- 
cies are induced by mechanical resonances in the deformable 
miiTor 

The coupling efficiency rj into a single-mode fiber is com- 
puted by taking into account two contributions: the tip-tilt er- 
ror which induces an offset of the stellar Airy pattern with 
respect to the fiber core (assuming the on-sky tip-tilt to be 
perfectly translated into an on-fiber jitter), and the Strehl ra- 
tio which represents the coherent part of the beam and can be 
taken as a multiplicative factor in th e computation of coupling 
efficiency ( Ruilier & Cassaing '2001'). The mean and RMS cou- 
pling efficiency after correction by adaptive optics are given in 
Table|2l 

The effect of coupling efficiency fluctuations, which pro- 
duce unequal intensities in the two arms of the interferome- 
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ter, is to induce a non-null transmission on the optical axis. 
Letting e/ be the relative intensity error between the two 
beams, the expression of the central transmission is T^{Q, 0) = 
ej/l6 (.OlhyierMigggi) . In order to keep the instrumental rejec- 
tion rate above the required 10^, relative intensity errors should 
be smaller than 1%, while the AO-corrected atmospheric tur- 
bulence induces intensity fluctuations as high as 8% in the L' 
band. A servo loop for intensity matching is therefore neces- 
sary to obtain a deep and stable null in the L' band. Such a 
device is not needed in the N band, where relative fluctuations 
of 1% are expected. 





Average 


Std deviation 


Tip-tilt [mas] 





14 


Strehl L' band 


0.80 


0.050 


Strehl N band 


0.97 


0.008 


Coupling L' band 


0.58 


0.044 


Coupling N band 


0.77 


0.008 



Table 2. RMS tip-tilt, mean and RMS Strehl ratio for an 
8-m telescope after connection by MACAO, as simulated by 
GENIEsim under typical atmospheric conditions at CeiTO 
Paranal (1" seeing, 11 m/s wind speed). The mean and RMS 
coupling efficiency into a single-mode fiber is deduced (taking 
into account the eff'ect of central obscuration). 



3.4. Scintillation 

Scintillation is the effect of rapid intensity fluctuations of a 
point-like source as a result of the interference of light rays 
diffracted by turbulent cells. The condition for scintillation to 
take place is that the Fresnel scale rp = (/I /; secz)'^^ is larger 
than the Fried scale ro, with h the height of the turbulent layer 
and z the zenith distance ( Qi^JIcilb'ich 1999). Even for high 
turbulent layers (20 km) and large zenith distances (z - 45°), 
weobtainrp ^ 0.3mwhilero - 1.1 mintheL' bandfor T'see- 
ing, showing that the Fried parameter is always larger than the 
Fresnel scale. Numerical estimates obtained with GENIEsim 
confirm this statement: scintillation does not induce intensity 
fluctuations larger than 0.2% in the L' band. Since ro cc A^^^ 
while rp <x A^^^, scintillation is also negligible for longer wave- 
lengths, including the N band. 

3.5. Polarization errors 

The polarization errors encountered by a nulling interferome- 
ter are of three main types: differential phase shift, differential 
attenuation and differential rotation can occur between the two 
beams for each linear polarization component. These effects 
can be translated into phase and intensity eiTors, but unlike 
phase and intensity errors due to turbulence, they are only due 
to the instrument itself. Polarization eiTors are thus expected 
to be mainly static or slowly drifting, even if instrumental vi- 
brations might induce some fluctuations. A highly symmetric 
design of the whole interferometer is required to reduce polar- 



ization issues as much as possible dSerabvn & Colavitdl200lh . 
In the following study, we will assume that the nulling instru- 
ment is designed so as to keep this contribution negligible in 
the instrumental leakage budget. This may not be the case in 
practice, because existing interferometers are not necessarily 
designed to meet the tight tolerances on polarization associated 
with nulling interferometry. This might force the nuller to op- 
erate in polarized light, by inserting polarizers at the beginning 
of the optical train^ and just before recombination. The instru- 
mental throughput would then be reduced by 50%, thereby in- 
creasing the required integration times. 

4. Real-time control and instrumental design 

We have seen that atmospheric turbulence causes additional 
stellar photons to leak through the null (instrumental leakage), 
and that this bias cannot be calibrated analytically as for ge- 
ometric leakage. Moreover, the stochastic fluctuation of in- 
strumental leakage is the source of an important noise con- 
tributi on, called systematic noise (iLav 2004) or variability 
noise JChazelas etalJl2005h . This section describes the three 
major control systems that are needed to reduce instrumental 
leakage and variability noise down to an acceptable level by 
correcting the effects of atmospheric turbulence in real time. As 
discussed in Sect.|3]and summarized in Table|3] an instrumen- 
tal rejection rate of 10^ requires the residual phase and intensity 
errors to be smaller than 4 nm and 1% RMS respectively. If this 
level of performance cannot be met, the use of calibrator stars 
will be further investigated to estimate and remove the contri- 
bution of instrumental leakage. 

4.1. Fringe tracking 

The purpose of a fringe tracker is to actively compensate for 
random phase fluctuations between the beams of the interfer- 
ometer at a given wavelength. A fringe sensing unit (FSU) 
measures the phase difference at a given sampling frequency 
and transmits this information to a controller which computes 
the compensation to be applied to one of the beams by an opti- 
cal delay line (ODL), as illustrated in Fig.|6] Both the FSU and 
the ODL introduce noise in the fringe tracking control loop, 
while the controller is (almost) noiseless. In order to be effi- 
cient, fringe tracking requires a high signal-to-noise ratio to 
be achieved on phase measurements on very short time scales 
(typically between 0.1 and 1 msec). It is therefore usually per- 
formed in the near-infrared (A < 2.4 fim) where stellar photons 
are numerous and the sky background rather low. 

Large interferometric facilities such as the VLTI or Keck 
Interferometer are (or will soon be) equipped with fringe track- 
ers. For instance, FINITO and PRIMA at VLTI will deliver 
a residual OPD of about 150 nm RMS on bright stars, re - 
spectively in the H and K bands dWilhelm & G itton"2003l). 
However, this is not sufficient to meet our goal on OPD stabil- 
ity, implying that a second stage of OPD control must be imple- 

^ In practice, a polarizer cannot be put in front of the primary or 
secondary mirror of a telescope, so that the first few mirrors of the 
optical train could still contribute to polarization errors. 
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Fig. 6. Simplified block-diagram of a fringe tracking system: 
the OPD error between two beams, characterized by its PSD 
■5 ATM, is sampled by a fringe sensing unit (FSU) feeding an 
OPD controller (C), which computes the correction to be ap- 
plied by the optical delay line (ODL). The latency represents 
the effect of the finite CPU and electronics time required 
between fringe sensing and correction, while the digital-to- 
analogue converter (DAC) represents the conversion from a 
digital command to an analogue OPD correction by the de- 
lay line. The OPD power spectrum after closed-loop control 
is denoted 5 res, while the noise introduced by the FSU and 
the ODL are respectively 5fsu and 5odl- 

mented inside the nulling instrument. In fact, the main limita- 
tion to the accuracy of the VLTI fringe tracker comes from the 
ODL, which does not efficiently correct for OPD fluctuations at 
frequencies higher than 50 Hz (see Fig.^J. The main require- 
ment for the additional fringe tracker is thus to be equipped 
with a low-noise ODL optimized for high repetition frequen- 
cies (up to 20 kHz). 



VLTI OPD control loop 
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Fig. 7. Simulation of the VLTI OPD control loop. The plot 
gives the PSDs at the input (before actuation: subscript "IN") 
and output (after actuation: subscript "OUT") of the loop for 
the three noise sources in the VLTI fringe tracking loop, i.e., 
atmospheric fluctuations (Atm, dotted curves), fringe sens- 
ing (FSU, dashed-dotted curves) and optical delay line (ODL, 
dashed curves). The total closed-loop OPD jitter (area under 
the solid line) is 133 nm RMS for a repetition frequency of 
4 kHz. 

We have simulated the performance of a fringe tracker op- 
timized to achieve small residual OPD on bright stars. The 



fringe tracker features a fringe sensor working both in the H 
and K band and uses a short-stroke fast piezo or voice-coil de- 
lay line to compensate the remaining OPD fluctuations. The 
control loop parameters are computed through a simultaneous 
optimization of the repetition frequency and of the controller 
parameters (a simple PID"*). The repetition frequency has been 
artificially limited to 20 kHz to comply with state-of-the-art 
back-end electronics. The delay line is assumed to have a per- 
fect response up to this frequency (short-stroke piezo transla- 
tors can achieve this kind of performance). The OPD control 
performance, summarized in Table|3 is computed in two cases: 

- the worst case corresponds to a -8/3 slope for the turbu- 
lence PSD at high frequencies (see Fig.O and includes the 
bimorph piston effect; 

- the best case takes into account the effect of pupil averaging 
on the PSD (-17/3 slope at high frequencies) and does not 
include the bimorph piston effect. 

The predicted performance ranges between 6 nm and 17 nm 
RMS . These figures could be improved by using a priori knowl- 
edge on the behaviour of atmospheric turbulence in the con- 
trol process, e.g. through Kalman filtering. First estimates have 
shown an improvement of about 2 nm on the residual OPD with 
respect to the above numbers (C. Petit, personal communica- 
tion). This performance, summarized in Table |3j is marginally 
compliant with the required 4 nm RMS. 

4.2. Compensation of longitudinal dispersion 

There are two different sources of dispersion to be corrected: 
the slowly varying differential phase due to unequal paths 
through wet air in the delay lines, and an additional rapid phase 
fluctuation due to column density fluctuations of dispersive 
components in the atmosphere. The first one is deterministic 
and easy to correct provided that the atmospheric conditions 
are monitored in the delay lines, while the second one is the 
real problem and is mainly associated with water vapour, as 
well as other molecular species such as ozone, carb on diox- 
ide, e tc, which are the main contributors below 2.5 jum llMathail 
l2004l) . 

In order to correct for both inter-band and intra-band dis- 
persion (Sect. 13. 2> . one should measure the differential col- 
umn density of water vapour experienced by the two beams 
in their paths through the atmosphere. This quantity can be in- 
ferred by measuring the phase difference between the beams at 
two different wavelengths, using the knowledge of the refrac- 
tive index of wet air across the infrared. In practice, a blind 
phase correction in the science band, relying for instance on H- 
and K-band measurements, is dangerous because it is not sensi- 
tive to possible biases, e.g. related to the imperfect knowledge 
of the refraction index or to the presence of other dispersive 
gases /M eisner & Le Pool3l2003h . Phase measurements in the 
science band are thus necessary to ensure a good co-phasing, 
using a part of the scientificafly useful signal. These measure- 
ments will be carried out at a slower rate, on one hand because 

* PID stands for "Proportional, Integral and Differential" and is a 
basic controller for closed-loop control. 
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of the lower signal-to-noise ratio on phase measurements at 
longer wavelengths and on the other hand because the possi- 
ble biases related to the extrapolation of the phase delay from 
the H/K bands towards longer wavelengths are mainly affected 
by global atmospheric parameters (temperature, partial pres- 
sures of gases, etc.) which are not expected to evolve quickly 
during the night. Another possibility would be to perform both 
fringe tracking and dispersion control by using science-band 
photons exclusively, in order to reduce the technical complexity 
of the instrument. This would however lead to slightly reduced 
performance, due to the worse signal-to-noise ratios on phase 
measurement at longer wavelengths, and would even become 
impractical in the N band. 

Interferometric facilities are currently not equipped with 
a dispersion compensation device, which should thus be inte- 
grated in the nuller itself. The required correction is two-fold: 
the inter-band phase error should be corrected with a dedicated 
short-stroke delay line (less than one fim RMS to be compen- 
sated), while intra-band dispersion could be corrected by intro- 
ducing a variable amount of dispersive material in the beams 
to cor rect for the slope of the phase inside the science wave- 
band d Koreskoetal .120031) . The performance of such a correc- 
tion system is summarized in Table|3lin the L' band, once again 
in two cases: 

- the worst case is obtained by assuming that phase measure- 
ments in the L' band are needed at each step of the control 
process, and assuming a -8/3 slope for the high-frequency 
asymptote of the atmospheric turbulence PSD. 

- the best case is computed assuming that H/K band phase 
measurements are sufficient for dispersion control and that 
the input PSD has a -17/3 slope at high frequencies. 

The closed-loop performance ranges between 4.4 nm and 
17 nm RMS for inter-band dispersion, between 1.0 nm and 
4.1 nm RMS for intra-band dispersion. It is marginally com- 
pliant with the requirement (4 nm RMS). 

4.3. Intensity matching 

The third real-time control loop to be implemented in the 
nulling instrument is a device to correct for unequal intensi- 
ties between the two beams. In fact, intensity fluctuations are 
not a problem per se: it is only the differential fluctuations be- 
tween the beams that matter. The most straightforward way to 
compensate for this effect is to measure the intensities of the 
two beams, determine which one is the brightest and reduce 
its intensity down to the level of the other one by means of 
an actuator This correction must take into account the effect 
of coupling into single-mode fibers, so that the intensity sen- 
sors must also be placed after the modal filters. Furthermore, it 
must be carried out in the science waveband, because coupling 
efficiency is wavelength-dependent. We will thus "waste" an 
additional part of the useful signal for intensity control. A pos- 
sible actuator could consist in a piezo-driven adjustable iris, 
which would reduce the beam size in a pupil plane while keep- 
ing it (almost) circular to preserve injection properties as much 
as possible. 
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The closed-loop control performance is summarized in 
Table |3]for the L' band, with a typical residual intensity mis- 
match of 4% when the loop is running at 1 kHz. This perfor- 
mance does not meet the required 1 % mismatch, which makes 
additional calibrations mandatory (see Sect. 15. 3> . Note that this 
control loop does not remove the fluctuations of the mean cou- 
pling efficiency (common to both beams), which has a typical 
standard deviation of 8%. 

4.4. Instrumental nulling performance 

Table |3] shows the nulling performance associated with the 
residual phase and intensity errors after closed-loop control 
for a G2V star at 20 pc. The RMS instrumental nulling ratio 
corresponds to the concept of variability noise ( Chaz elas et alJ 
i2005i) . Its power spectrum is illustrated in Fig. |S] showing 
an almost flat spectrum for frequencies below 1 kHz and a 
rapid decrease of the power content above 1 kHz. Variability 
noise thus behaves as a white noise at typical read-out fre- 
quencies (~10 Hz), so that its cumulated contribution in- 
creases as the square root of integration time just as for shot 
noise. We do not take into ac count here the possible 1 //-type 
noises dChazelas et al.lEoo3) . which might alter the slope of 
the PSD at low frequencies and thus change the behaviour of 
long-term variability noise by introducing more power at low 
frequencies. This assumption is supported by the use of real- 
time control loops, which will help remove possible drifts in 
the instrumental response. 





worst case 


best case 


goal 


Piston 


17nm @ 20kHz 


6.2nm @ 13kHz 


< 4nm 


Inter-band 


17nm @ 200Hz 


4.4nm @ 300Hz 


< 4nm 


Intra-band 


4.1nm @ 200Hz 


l.Onm @ 300Hz 


< 4nm 


Intensity 


4% @ IkHz 


4% @ IkHz 


< 1% 


Total null 


9.7 X 10-" 


6.2 X 10-* 


4.8 X 10-" 


Instr. null 


5.0 X 10-* 


1.5 X lO-'' 


10-5 


RMS null 


5.8 X 10-* 


1.6 X 10-" 


10-5 



Table 3. Performance of an L' band nuller at VLTI as simu- 
lated with GENIEsim on a 100 sec observation block, taking 
into account all turbulence-induced errors. The performance of 
the control loops are summarized together with their repetition 
frequencies in pessimistic and optimistic cases (see text). The 
target performance discussed in Sect.|3appears in the last col- 
umn. The total null is the mean nulling ratio including both 
the geometric and instrumental leakage contributions, and is 
computed as the ratio between the total input and output stellar 
fluxes (from both telescopes). The last line gives the standard 
deviation of the instrumental nulling ratio for this 100 sec ob- 
servation block. 



5. Post-processing of nulling data 

Because the expected performance of state-of-the-art real-time 
control loops is not sufficient to reduce instrumental leakage 
down to the required level (10-^ in the L' band), and because 
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Fig. 8. Power spectrum of instrumental nulling due to phase 
and intensity eiTors after closed-loop control, obtained for a 
100 sec simulation with GENlEsim in the L' band. The semi- 
analytical si mulat i on tak es into account the bilinear error terms 
discussed by Lav (2004), but not the possible 1 //-type noises 
discussed by Chazelas et al. (2005.) . The power spectrum has 
been smoothed with a standard IDL routine to show its mean 
behaviour: almost flat at low frequencies, while the logarithmic 
slope at high frequencies is about -2. 

other sources such as geometric leakage and background also 
contribute to the detected signal after destructive beam combi- 
nation, the exozodiacal disk signal can only be extracted from 
the raw data after a final and critical phase of coiTection: post- 
processing. 

5.1. Background subtraction 

The first step in the data analysis procedure is to remove the 
contribution of the background. In order to extract the use- 
ful data from the incoherent background in the nulled data, 
modulation techniques are required. The classical "chopping- 
nodding" technique consists in tilting the secondary mirror by 
a small quantity (a few arcseconds) in order to measure the 
background in an empty sky region close to the source and 
subtract it from the on-source measurement. The chopping se- 
quence is repeated at a few Hertz. Nodding consists in repeat- 
ing the whole chopping sequence off'-source, by slightly tilting 
the telescope, in order to remove the possible gradient in the 
background emission. 

Because of the large temporal fluctuations of the back- 
ground emission, chopping at a few Hz is not sufficient to re- 
duce sky noise below the level of the scientifically useful sig- 
nal after background subtraction. Other techniques are there- 
fore considered in order to measure the background emission 
simultaneously to the observations, e.g. by inserting additional 
fibers in the transmitted field-of-view of the interferometer. 
At VLTI for instance, the use of variable curvature mirrors in 
the delay lines provides a clean field-of-view 2" in diameter, 
which allows us to reduce the contribution of stellar light in 
the background fibers down to about 10 of the background 
contribution. Another solution for background subtraction is to 
use phase modulation techniques, which require the telescope 



pupils to be divided into two parts in order to modula te be- 
tween two different nulled outputs JSerabvn et alJl2004 . This 
is the only practical solution for an N-band nuller, for which the 
contrast between background and exozodiacal fluxes is huge 
(almost 10^ see TablegJ. 

5.2. Geometric ieakage calibration 

An important step in the post-processing of nulling interfer- 
ometry data is the calibration of geometric leakage. We have 
discussed in Sect. l2.6l that the residual contribution of geomet- 
ric leakage after calibration should not exceed 10"^ of the ini- 
tial stellar flux in the L' band, while typical rejection rates of 
2000 are expected (TableQ. The calibration of geometric leak- 
age, based on the analytical expression (|3} of the rejection rate, 
should therefore reach a precision of about 2% in order to re- 
trieve the useful exozodiacal disk signal with a good signal-to- 
noise ratio. Closer stars, for which the rejection rate is smaller, 
require even better calibration accuracies. A typical precision 
of 1% for geometric leakage calibration is thus required in the 
U band, while such a calibration is marginally needed in the 
N band since the achievable rejection rate (~16000) already 
matches the requirement on starlight rejection (10 

Assuming that the baseline B and the effective wavelength 
A are known with a good accuracy, a good knowledge of 
the stellar angular radius 0, and limb-darkening coefficient u^ 
is sufficient in order to compute the geometric leakage. The 
global optical transmission of the interferometer must also be 
known in order to convert this value into a number of photo- 
electrons measured at the detector While the interferometric 
transmission can be straightforwardly measured with a high 
precision at the constructive output of the nuller where most 
of the stellar photons end up, the imperfect knowledge of stel- 
lar diameters has a large (and generally dominant) contribution 
to the caUbration error budget. Differentiating equation (|3j, we 
obtain a required precision A0,/6t = 0.5AN/N - 0.5% on the 
stellar diameter knowledge. 

The most precise way to determine stellar angular radii is 
currently based on stellar interferometry. Stellar diameter mea- 
surements with an accuracy as good as 0.2% have bee n demon- 
strate d on resolved stars such as a Cen A and B ( Kervella et alJ 
l2003h . However, the Darwin/TPF-1 targets, with typical an- 
gular diameters of 1 mas and below, are only marginally re- 
solved even with the longest baselines cuiTently available: the 
AMBER instrument at VLTI provides at best an angular reso- 
lution of 1 .6 mas in J band. A precision of a few 0. 1 % on diam- 
eter measurements with AMBER could thus only be achieved 
for the closest stars in the Darwin catalogue (< 10 pc). An 
additional difficulty comes from the extrapolation of the limb- 
darkening law from the J band to the L' band, since limb- 
darkeni n^ coefficient s are currently modelled only up to the 
K band (Claret 2000). Lunar occultation measurements could 
also be used to determine the diameters of the brightest targets 
of the Darwin catalogue {K ~ 3) with a limiting resolution of 
~ 1 mas, but with the restriction that only about 10% of the sky 
faUs on the Moon's path (Fors et al.-2004.) . 
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For stars more distant than about 10 pc, for which inter- 
ferometric diameter measurements will probably not reach a 
0.5% precision, indirect measurements based on surface bright- 
ness relations might be used. Empirical laws make it possi- 
ble to predict the limb-darkened angular diameters of dwarfs 
and subgiants using their d ereddened Johnson magnitudes, or 
their effective temperature aKervella et aljE004 ). Intrinsic dis- 
persions as small as 0.3% are obtained for the relations based 
on B-L magnitudes. However, the apparent magnitudes of the 
Darwin/TPF-1 targets need to be accurately known (to better 
than 0.02 mag): this is usually not the case at infrared colours 
and might mandate a dedicated photometry programme. 

All in all, it is expected that diameter measurements with a 
precision of 0.5% will be obtained on most of the Darwin/TPF- 
I targets. A 1% precision is still considered as a worst case 
scenario in the performance budget. 

5.3. Instrumental leakage calibration 

After calibration of geometric leakage, there is one remaining 
step of calibration to be carried out: the evaluation and sub- 
traction of the mean instrumental leakage from the nulled data. 
Calibration of the instrumental response is a routine procedure 
in stellar interferometry, and is achieved by observing a well- 
know calibra tor star, usually c hosen as a late-type unresolved 
giant star (M erand et aklbOOSl) . These measurements are car- 
ried out just before and after the observation of the target star, 
and as close as possible to it in the sky, so that the atmospheric 
properties do not change significantly. 

This method can be extended to measure instrumental leak- 
age in nulling interferometry. In this case, calibrator stars must 
be carefully chosen: even if instrumental leakage does not de- 
pend on the angular diameter of the star to the first order, the 
control loop performance strongly depends on the stellar flux 
and spectrum, so that the calibrator star should in fact have the 
same flux and spectrum as the target star in order to get the 
same instrumental leakage in both cases. Having the same flux 
and spectrum means that the angular diameters should also be 
the same, and consequently that the calibrator star will pro- 
duce the same amount of geometric leakage as the target star. 
A precise knowledge of the calibrator's angular diameter is thus 
needed in order to remove the contribution of geometric leak- 
age and obtain a good estimation of instrumental leakage. 

The final estimation Z{A) of the exozodiacal signal after 
both geometric and instrumental leakage calibrations can be 
expressed as follows: 

Z{A) = S,{A) - -j^F,{A) - \S,{A) - -j^F,{A)\ , (5) 

with Si{A) and 5c(/l) the background-subtracted nulled fluxes 
for the target and calibrator stars, 0, and 6c their angular radii, 
Fi{A) and f ^(/l) their fluxes, measured at the constructive out- 
put of the interferometer. The first subtracted term represents 
the contribution of geometric leakage while the second one 
(between brackets) is the estimation of instrumental leakage 
obtained with the calibrator star. 



6. Detection of exozodiacal dust disks witli GENIE 

In the following paragraphs, we evaluate the sensitivity of 
GENIE, the Ground-based European Nulling Interferometer 
Experiment jGondoin et al J2004h . to exozodiacal dust for typ- 
ical Darwin/TPF-I targets. 

6.1. Background-limited signal-to-noise ratio 

Before computing the actual detection performance of the 
nuller, let us assess the required integration times in a non- 
turbulent environment, where the only source of noise is the 
shot noise from all sources inside the field-of-view. The pho- 
ton budget of Table 0]includes the geometric and instrumental 
leakage contributions given in Table |3] (worst case scenario), 
the contribution of the dust disk itself^ and the diff'erent sources 
of background emission (sky and instrument). 

We assume an overall instrumental throughput (including 
VLTI, GENIE and detectors) of 1 % both in the U and N bands. 
Table 0]gives the time required to detect a 20-zodi disk around 
a G2V star at 20 pc with a signal-to-noise ratio of 5. While the 
integration time in the L' band is reasonable (15 minutes), it is 
very large in the N band (11 hours). In fact, the difficulty in the 
N band comes not only from the required integration time, but 
also from the huge contrast between the exozodiacal and back- 
ground fluxes (almost 10^) at the output of the nuller These two 
background-related issues are the main reasons for the choic e 
of the U band in the context of GENIE JGondoin et alJl2004 . 
Note that, even in the L' band, the required background calibra- 
tion accuracy is challenging: in order to reach a final signal-to- 
noise ratio larger than 5, the background must be removed with 
an accuracy of 10"^, which is within reach using the advanced 
background subtraction methods discussed in Sect. 15.1! 

6.2. Final signal-to-noise ratio 

In this section, we evaluate the global sensitivity of the GENIE 
instrument by simulating a 30-minute observation of a Sun-like 
star at 20 pc surrounded by a 20-zodi disk. In Table |5l we first 
give the raw sensitivity of the nuller, i.e., without calibration 
of instrumental leakage (very conservative approach assuming 
that no convenient calibrator star can be found). In this case, 
the mean instrumental leakage fully contributes as a bias which 
adds to the exozodiacal disk signal that we want to detect. This 
implies that the mean instrumental leakage should be at least 
five times smaller than the disk contribution to ensure a safe 
detection. Table |5] shows that, depending on the model for at- 
mospheric turbulence, the raw sensitivity ranges between 530 
and 1800 zodis for a signal-to-noise ratio of 5 across the whole 
L' band. Instrumental leakage is the main contributor to the 
noise budget in that case. Note that shot noise, detector noise 
and variability noise are almost negligible after only 30 min of 
integration. This would not be the case if spectral dispersion 
was used: in order to keep the same sensitivity, the integration 

^ The part of the exozodiacal emission transmitted by the interfer- 
ence fringe pattern shown in Fig.Qtypically ranges between 30% and 
40%. 
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L' 


N 


Total stpllflr 'iicTifil FTvl 


3.1 


0.51 


Total disk signal [Jy] 


3.0 X 10""' 


5.1 X 10"^ 


Sky brightness [Jy/as"] 


5.0 


690 


VLTI brightness [Jy/as^] 


115 


21300 


GENIE brightness [Jy/as-] 


29 


5320 


Total bckg signal [Jy] 


1.8 


2550 


Final stellar leakage [el/s] 


1.0 X lO'' 


1.5 X 10^ 


Final disk signal [el/s] 


3.6 X 10^ 


4.1 X 10^ 


Final bckg signal [el/s] 


4.8 X 10* 


2.7 X 10'" 


Shot Noise [el/s"-] 


2.2 X 10^ 


1.7 X 10^ 


Time for SNR = 5 [sec] 


950 


40000 



Table 4. In the first part of the table, we compare the star and 
exozodiacal disk fluxes to the thermal background, computed 
from its individual contributors (Paranal night sky, VLTI opti- 
cal train and GENIE instrument) by integration over the field- 
of-view defined by a single-mode fiber (Q = A^/S). The target 
is a Sun-like star at 20 pc surrounded by a 20-zodi disk. The 
values are given in Jy at the central wavelength of the infrared 
atmospheric windows and "at the entrance of the Earth's at- 
mosphere" (an equivalent emission is given for the background 
contributors as if they originated from outside the Earth atmo- 
sphere). Second part: signals in electrons per second actually 
detected at the GENIE science detector. The total shot noise is 
deduced, and used to compute the time needed for a photomet- 
ric signal-to-noise ratio of 5. 



time should then be increased by a factor equal to the square 
root of the number of spectral channels. 

In a second step, we compute the performance assuming 
that an appropriate calibrator star can be found, i.e., with simi- 
lar flux, spectrum and angular diameter than the target star but 
without circumsteUar dust. In this case, the calibration preci- 
sion for instrumental leakage is mainly limited by the imper- 
fect knowledge of the calibrator angular diameter, assumed to 
be the same than for the target star. This contribution is dom- 
inant in the noise budget, so that the GENIE sensitivity could 
benefit from improvements in the precision of stellar diame- 
ter measurements. In practice, the calibrator star will not be 
perfectly identical to the target, inducing a bias in the calibra- 
tion process due to the different behaviour of the control loops 
for different stellar magnitude and/or spectra. This bias is not 
dominant in the calibration process as long as the H, K and L' 
magnitudes of the target and calibrator do not differ by more 
than about 0.2 magnitudes. A convenient calibrator will proba- 
bly not be available in all cases, so that the actual performance 
of the nuller could be somewhat degraded with respect to the 
final values of Table|5](Absil et al., in preparation). 

6.3. Influence of stellar type and distance 

In order to assess the scientific pertinence of GENIE in the con- 
text of the Darwin/TPF-I programmes, we have evaluated the 
detection performance for exozodiacal dust clouds around typ- 
ical Darwin/TPF-I targets, i.e. late-type dwarfs with distances 
ranging between 5 and 25 pc. In Fig.|9lare plotted the expected 





worst case 


best case 


9n-7nHi sionfil fp-l 


6.5 X lO*" 


8.7 X 10' 


Shot noise [e-] 


1.3 X 10' 


1.5 X 10' 


Detector noise [e-] 


3.4 X 10^ 


4.0 X 10^ 


Variablility noise [e-] 


1.9 X lO'* 


5.9 X 10'' 


Calibrated geom. leakage [e-] 


2.2 X 10' 


1.5 X 10' 


Raw instr. leakage [e-] 


1.2 X 10' 


4.6 X 10*^ 


Calibrated instr. leakage [e-] 


3.8 X 10' 


2.1 X 10' 


Zodis for SNR=5 (raw) 


1800 


530 


Zodis for SNR=5 (calibrated) 


56 


34 



Table 5. Expected sensitivity of the GENIE instrument at 
VLTI, given in number of zodis that can be detected around 
a Sun-like at 20 pc in 30 min in the full L' band using the UT2- 
UT3 baseline (47 m). The worst case corresponds to conser- 
vative control loop performance and 1% precision on diame- 
ter measurements, while the best case relies on optimistic con- 
trol loop performance and 0.5% precision on angular diame- 
ters. Each individual contribution is given in photo-electrons 
detected at the nulled output. The raw sensitivity is given with- 
out any calibration of instrumental leakage, while the calibrated 
sensitivity assumes that a calibrator star with similar character- 
istics as the target star is used to perform instrumental leakage 
calibration. 



detection levels for the GENIE instrument with 30 minutes of 
integration on the 8-m UTs, while Fig. ^| gives the perfor- 
mance on the 1.8-m ATs for an integration time of 2 hours. 
The baseline lengths have been optimized in order to reach an 
optimum sensitivity to exozodiacal dust, within the limitations 
imposed by the VLT interferometer (UT-UT baselines com- 
prised between 47 and 130 m, AT- AT baselines ranging from 
8 to 202 m). Here again, we have assumed that a calibrator 
star similar to the target star can be found for the evaluation 
of instrumental leakage. This optimistic assumption is at least 
partially compensated by the fact that we use conservative as- 
sumptions on atmospheric turbulence properties and on stellar 
diameter knowledge (cf. the "worst case" of Table|5}. 

The performance estimates in the case of UTs show a maxi- 
mum efliciency for exozodiacal dust detection at 25 pc for FOV, 
20 pc for GOV, 15 pc for KOV and 10 pc for IVIOV stars. This 
maximum is reached when the two main sources of noise are 
well balanced: for nearby targets, which are partially resolved, 
the contribution of geometric leakage is dominant, while for 
more distant and thus fainter targets the main contribution 
comes from shot noise. The strong limitation on the sensitivity 
for close targets can be overcome by using smaller baselines, 
which are available only for AT pairs at the VLTI. This is illus- 
trated in Fig. [TO] with the shortest baselines, one can achieve 
exozodiacal dust detection down to the 20-zodi level for bright 
and close targets (especially F-type stars), which are however 
scarce in the Darwin star catalogue. 

These simulations indicate that exozodi detection down to 
a density level of about 50 times our local zodiacal cloud is a 
realistic goal for most F and G stars in the Darwin catalogue, 
while exozodiacal densities of about 100 zodis could be de- 
tected around K-type stars. In the case of M stars, only the 
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closer stars could lead to a pertinent detection level, of about 
200 zodis. 
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Fig. 9. Sensitivity of GENIE to dust disks around nearby main- 
sequence stars, expressed in units relative to the Solar zodiacal 
cloud. The simulations have been obtained with optimized UT- 
UT baselines for 30 minute exposures. The smallest baseline 
(47 m on UT2-UT3) is used in most cases because the dominant 
noise source generally comes from the calibration of geometric 
leakage for both the target and calibrator stars. 
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Fig. 10. Sensitivity of GENIE to dust disks around nearby 
main-sequence stars, expressed in units relative to the Solar 
zodiacal cloud. The simulations have been obtained with op- 
timized AT-AT baselines in 2 hours of integration. The avail- 
ability of very small baselines (down to 8 m) allows us to re- 
duce the contribution of geometric leakage and thus improves 
the detection levels of Fig.|51at short distances. 



7. Conclusions and perspectives 

In this paper, we have demonstrated the exceptional capabili- 
ties of ground-based nulling interferometers to detect faint cir- 
cumstellar features such as exozodiacal disks. Ground-based 
nulling interferometry requires high-performance closed-loop 
control of atmospheric turbulence in the L' band, and ex- 
tremely precise background subtraction in the N band. In this 



paper, we have computed the expected performance of state- 
of-the-art servo loops for piston, dispersion and intensity con- 
trol, allowing the stabilization of the instrumental nulling ra- 
tio. Calibration procedures have been discussed to further re- 
move stellar light from the nulled output of the instrument. 
This comprehensive study shows that an L'-band nuller such 
as GENIE at VLTI could detect exozodiacal disks about 50 to 
100 times as dense as the Solar zodiacal cloud. This instru- 
ment will thus significantly improve the detection performance 
of current infrared and sub-millimetric facilities, and above all 
will peer into the cuiTently inaccessible warm region of dust 
disks located within a few AU of nearby late-type dwarfs (typ- 
ical Darwin/TPF-I targets). 

The simulated performance of the nulling instrument criti- 
cally depends on a number of noise contributions which could 
be improved in different ways. Geometrical leakage could be 
reduced by using short baselines, typically ranging between 8 
and 40 m, while its calibration would be improved by obtain- 
ing accurate stellar diameter measurements at long baselines 
(> 300 m) and/or precise photometry to derive diameters from 
surface-brightness relationships, and by extending the limb- 
darkening models towards the mid-IR wavelengths used by the 
nulling instruments. Instrumental leakage depends strongly on 
the ability to correct the effects of atmospheric turbulence, and 
would therefore benefit from advances in real-time control al- 
gorithms and compensation devices. Polarization errors, which 
have been discussed only briefly in this paper, could also be 
critical and need to be carefully addressed at the design level. 
Finally, advanced background subtraction methods such as the 
phase modulations techniques foreseen for the Darwin/TPF-I 
missions might further improve the overall performance of a 
ground-based nuller. 

Acknowledgements. The authors wish to thank their ESA and ESO 
colleagues, the two industrial partners (Alcatel Space and EADS- 
Astrium) and the GENIE scientific team for their support and con- 
tributions. The authors are also grateful to the anonymous referee and 
to Jean Surdej for their precious help in improving the quality and 
readability of the paper. 



References 

Absil, O. 2001, Diploma thesis. University of Liege, 
http://vela.astro.ulg.ac.be/themes/telins/harigs/ 

Absil, O., Bakker, E., Scholler, M., & Gondoin, P 2004, 
in Proc. SPIE, Vol. 5491, New Frontiers in Stellar 
Interferometry, ed. W. Traub, 1320 

Absil, O., den Hartog, R., Erd, C, et al. 2003a, in Toward Other 
Earths: Darwin/TPF and the Search for Extrasolar Terrestrial 
Planets, Vol. SP-539 (ESA's PubUcation Division), 317 

Absil, O., Kaltenegger, L., Eiroa, C, et al. 2003b, 
in Toward Other Earths: Darwin/TPF and the Search 
for Exti-asolar Terrestrial Planets, Vol. SP-539 (ESA's 
Publication Division), 323 

Absil, O., Karlsson, A., & Kaltenegger, L. 2003c, in Proc. 
SPIE, Vol. 4852, Interferometry in Space, ed. M. Shao, 431 

Allen, D. & Barton, J. 1981, PASP, 93, 381 



14 



Absil et al.: Performance study of ground-based infrared Bracewell interferometers 



Angel, J. 1989, in The Next Generation Space Telescope, 

ed. P. Bely, C. Burrows, & G. IlUngworth (HST Sc. I., 

Baltimore), 81 
Angel, J. & Woolf, N. 1997, ApJ, 475, 373 
Arsenault, R., Alonso, J., Bonnet, H., et al. 2003, in Proc. 

SPIE, Vol. 4839, Adaptive Optical System Technologies, ed. 

P. Wizinowich & D. Bonaccini, 174 
Beichman, C, Bryden, G., Rieke, G., et al. 2005, ApJ, 622, 

1160 

Beichman, C, Woolf, N., & Lindensmith, C. 1999, TPF: A 
NASA Origins Program to Search for Habitable Planets (Jet 
Propulsion Laboratory: Pasadena, CA), JPL 99-3 

Bracewell, R. 1978, Nature, 274, 780 

Chazelas, B., Brachet, R, Borde, P., et al. 2005, Applied Optics, 

accepted 
Claret, A. 2000, A&A, 363, 1081 

Claret, A., Diaz Cordoves, J., & Gimenez, A. 1995, A&AS, 
1 14, 247 

Colavita, M., Wizinovich, P., & Akeson, R. 2004, in Proc. 

SPIE, Vol. 5491, New Frontiers in Stellar Interferometry, ed. 

W. Traub, 454-463 
Conan, J.-M., Rousset, G., & Madec, P-Y. 1995, J. Opt. Soc. 

Am. A, 12, 1559 
Dermott, S., Grogan, K., Holmes, E., & Wyatt, M. 1998, 

in Proceedings of the Exo-zodiacal Dust Workshop, ed. 

D. Backman, L. CarofF, S. Sanford, & D. Wooden, Vol. 1998- 

10155 (NASA/CP), 59-63 
Fors, O., Richichi, A., Nunez, J., & Prades, A. 2004, A&A, 

419, 285 

Fridlund, M. 2000, ESA bulletin, 103, 20 

Ghndemann, A., Albertsen, M., Andolfato, L., et al. 2004, 

in Proc. SPIE, Vol. 5491, New Frontiers in Stellar 

Interferometry, ed. W. Traub, 447^53 
Gondoin, P., Absil, O., den Hartog, R., et al. 2004, in Proc. 

SPIE, Vol. 5491, New Frontiers in Stellar Interferometry, ed. 

W. Traub, 775-786 
Herbst, T. & Hinz, P 2004, in Proc. SPIE, Vol. 5491, New 

Frontiers in Stellar Interferometry, ed. W. Traub, 383-390 
Karisson, A. & Mennesson, B. 2000, in Proc. SPIE, Vol. 

4006, Interferometry in Optical Astronomy, ed. P. Lena & 

A. Quirrenbach, 871 
Kelsall, T., Weiland, J., Franz, B., et al. 1998, ApJ, 508, 44 
Kervella, P., Thevenin, F., di Folco, E., & Segransan, D. 2004, 

A&A, 426, 297 

Kervella, P, Thevenin, F., Segransan, D., et al. 2003, A&A, 
404, 1087 

Kolmogorov, A. 1941, Doklady Akad. Nauk. S.S.S.R., 30, 229 
Koresko, C, Mennesson, B., Serbyn, E., et al. 2003, in Proc 

SPIE, Vol. 4838, Interferometry for Optical Astronomy II, 

ed. W. Traub, 625 
Kuchner, M. & Brown, M. 2000, PASP, 1 12, 827 
Kuchner,M., Brown, M., & Koresko, C. 1998, PASP, 1 10, 1336 
Lawson, P. & Dooley, J. 2005, Technology Plan for the 

Terrestrial Planet Finder Interferometer (Jet Propulsion 

Laboratory: Pasadena, CA), JPL 05-5 
Lay, O. 1997, A&AS, 122, 535 
Lay, O. 2004, AppUed Optics, 43, 6100 



Leger, A., Mariotti, J., Mennesson, B., et al. 1996, Icarus, 123, 
149 

Mathar, R. 2004, Applied Optics, 43, 928 

Meisner, J. & Le Poole, R. 2003, in Proc. SPIE, Vol. 4838, 

Interferometry in Optical Astronomy n, ed. W. Traub, 609 
Mennesson, B., OUivier, M., & RuiUer, C. 2002, J. Opt. Soc. 

Am. A, 19, 596 

Merand, A., Borde, P, & Coude du Foresto, V. 2005, A&A, 

433, 1155 

Ollivier, M. 1999, PhD thesis, Universite de Paris XI 
Quirrenbach, A. 1999, in Principles of Long BaseUne Stellar 

Interferometry, ed. P. Lawson (JPL publ.), 71 
Roddier, F. 1981, in Progress in optics, ed. E. Wolf, Vol. XIX 

(North-Holland, Amsterdam), 281 
Ruilier, C. & Cassaing, F 2001, J. Opt. Soc. Am. A., 18, 143 
Serabyn, E., Booth, A., Colavita, M., et al. 2004, in Proc. 

SPIE, Vol. 5491, New Frontiers in Stellar Interferometry, ed. 

W. Traub, 806 

Serabyn, E. & Colavita, M. 2001, Apphed Optics, 40, 1668 
Stankov, A., Kaltenegger, L., & Eiroa, C. 2005, ESA/ESTEC 

internal publication DMS/SCI-A/DARWIN/289 
ten Brummelaar, T. 1995, Optics Communications, 1 15, 417 
Verinaud, C. & Cassaing, F. 2001, A&A, 365, 314 
Wilhelm, R. & Gitton, P 2003, in Toward Other Earths: 
Darwin/TPF and the Search for Extrasolar Terrestrial 
Planets, Vol. SP-539 (ESA's Publication Division), 659 
Young, A., Milone, E., & Stagg, C. 1994, A&AS, 105, 259 



